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SUMMARY 


The problem of eigenvalue assignment in a linear time-invariant system 
using output feedback is considered. New sufficient conditions are derived to 
assign an almost arbitrary set of minimum (n,m + r - 1) distinct eigenvalues 
where n, m, and r are the number of states, inputs, and outputs, respec- 
tively. These conditions precisely identify. the class of systems where such an 
assignment is impossible. The synthesis technique also highlights the freedom 
in selection of closed-loop eigenvectors under output feedback. The utility of 
eigenvalue/eigenvector assignment in transient response shaping is illustrated 
by the design of a controller for the lateral dynamics of an aircraft. 


INTRODUCTION 

Control system design based on eigenvalue or pole assignment has received 
a great deal of attention in the literature. It is well known that for a con- 
trollable system, if state variable feedback is employed, all eigenvalues can be 
assigned (ref. 1). Also it is known that for multi-input systems, the feedback 
law assigning a given set of eigenvalues is not unique and that different con- 
trol laws can yield identical eigenvalues while yielding radically different 
eigenvectors. Since the eigenvectors determine the influence of each eigenvalue 
on each state variable response, failure to use the multi-input design freedom 
fully may result in undesirable mode coupling and other poor transient behavior. 
For n-state feedback systems, it has been shown (ref. 2) that with m inputs, 
in addition to the assignment of all n eigenvalues, up to m entries in each 
eigenvector can be arbitrarily assigned. However, the problem of eigenvalue 
assignment using output feedback instead of state feedback has not yet been com- 
pletely resolved. The problem of determining conditions under which all eigen- 
values of a system can be arbitrarily assigned to a system under output feedback 
has been investigated in references 3 and 4. Bounds on the number of states, in 
terms of number of inputs, outputs, controllability, and observability indices 
are established for complete pole assignability. Reference 5 shows that for a 
system with r outputs, if mr > n then the system is pole-assignable provided 
the feedback gain elements are allowed to be complex numbers. References 6 and 
7 address the converse problem: given a controllable, observable system, how 

many eigenvalues can be arbitrarily assigned to the system. In general, it is 
concluded (ref. 7) that minimum (n,m + r - 1) eigenvalues can "almost" always 
be assigned to the system using output feedback. The qualification "almost" was 
introduced to cover classes of systems where such an assignment is impossible. 

In effect, the analysis in reference 7 does not precisely determine the condi- 
tions under which (m + r - 1) eigenvalues cannot be assigned to the system. 

This report considers the problem of determining the number of eigenvalues 
assignable to a given system. By formulating an eigenvalue/eigenvector assign- 
ment problem, sufficient conditions required for the assignment of minimum 
(n,m + r - 1) eigenvalues are derived. These conditions precisely identify 


the class of systems which can be assigned only d eigenvalues, where maximum 
(m,r) <d < (m+r - 1). The new formulation permits the development of an 
algorithm to assign (m + r - 1) eigenvalues. In addition, (r - 1) eigenvectors 
can be partially assigned with, at most, m entries in each vector arbitrarily 
chosen. In the event n>(m+r-1), various synthesis alternatives to stabi- 
lize the system are also investigated since in this case all system eigenvalues 
cannot be assigned. The counter example of reference 7 is used to demonstrate 
the utility of the new sufficient conditions in identifying systems which cannot 
be assigned (m + r - 1) eigenvalues. Finally, the advantage of both eigenvalue 
and eigenvector assignments in response shaping is illustrated by designing a 
controller to meet the lateral handling qualities specifications for an aircraft 


SYMBOLS 

Values are given in SI and U.S. Customary Units. Calculations were made in 
U.S. Customary Units. 

A system matrix 

A € Rn>fn n X n real matrix A 

3y lateral acceleration, m/sec2 (ft/sec^) 

B input matrix 

C measurement matrix 

c measurement vector 

5 matrix used in equation (45) 

matrices defined by equation (B4) 
d,i,j,k,t indices 

vector defined in notation ( 1 ) in appendix B 
matrix defined by equation (9) 
vector defined by step 3(a) of appendix B 
matrix defined by equation (10) 
vector defined by step 3(a) of appendix B 

matrices used in equation (23) and derived by substituting equa- 
tion (6) into equation (7) 

vector defined by step 3(a) of appendix B 

( index) order identity matrix 

2 


fj(k-l) 


Sj 


(k) 




(k-1) 


^index 



j e (A(k)} 


j is an element of index set 


K feedback matrix 

matrices used in equation (23) and derived by substituting equa- 
tion (6) into equation (7) 

I feedback vector defined by equation (17) 

mO^) ^ ^ ) matrices defined by notation (6) in appendix B 

m number of inputs 

max maximum value 

min minimum value 

N matrix defined by equation (13) 

n number of states 

P matrix used in equation (46) and defined immediately afterward 

p roll rate, deg/sec 

Q(k-1)^Q(i) matrices defined by equation (B3) 

q( 0) matrix defined by notation (3) in appendix B 

q feedback vector defined by equation (18) 

r number of outputs 

S matrix defined by equation (8) 

s vector defined immediately prior to equation (20) 

To,T'| transformation matrices defined where used 

u input vector defined by equation (2) 

V modal matrix (matrix of eigenvectors) 

v(n) matrix used in equation (12) 

V ,Vj^ eigenvectors 

w vector used in equation (6); a partition of v 

X state vector defined by equations (1) 

X G Rn n X 1 real vector x 
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y output vector defined by equations (1) 

matrices containing the vectors z^; i = 1 to r and 1 to t, 
respectively 

z vector used in equation (6); a partition of v 

6 sideslip, deg 

aO ) , Aj. ( 1 1 ) , set of indices defined in appendix B 

aileron angular deflection, deg 
determinant defined in step 1 of appendix B 
rudder angular deflection, deg 

eXjj,ezij perturbed quantities defined in step 1 and step 4 of appendix B, 
respectively 

0 vector defined by notation (6) in appendix B 

Ap,A^. diagonal eigenvalue matrices used in equations (12) and (B6), 
respectively 

X,Xi,Xj^ eigenvalues 

^,ri,Pl,P2 vectors defined by equations (21) 

djj scalar defined by notation (4) in appendix B 

(|) bank angle, deg 

i> yaw rate, deg/sec 

Superscripts: 

-1 matrix inverse 

transformed quantity 

Upper case letters of the alphabet indicate matrices; matrix subscripts 
indicate partitioned quantities. Dot over a quantity denotes derivative with 
respect to time. Prime denotes transpose. 

EIGENVALUE/EIGENVECTOR ASSIGNMENT FORMULATION 

Consider a linear, time invariant, multivariable, controllable, observable 
system 
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X = Ax + Bu 


(1) 


y = Cx 


where x ^ u e rhIj y e pr^ and B and C are full rank; and for a non- 

trivial problem formulation, assume m,r > 1 and m,r < n. The problem is to 
find a control law of the form 

u = Ky (2) 

in order to assign arbitrary eigenvalues for the closed-loop system. To indi- 
cate clearly the freedom available in the selection of closed-loop eigenvalues 
and eigenvectors under output feedback, the measurement matrix C is assumed 
to be in a special canonical form: 


C = [Cl : C2] 


Cl = 



C2 = 


0 

c 




(3) 


c = [1 1 . . . 1]J 

Here, t = r - 1, C-| s c e Rl><n-t^ denotes a t'^*’ order identity 

matrix. Appendix A details a procedure for reducing any system (C,A,B) to this 
special form. 


The closed-loop system matrix (A + BKC) after applying feedback law (2) 
satisfies 


(A + BKOvi = X^Vi (i = 1 , 2, . . . , n) (4) 

where Xj^ is the ith eigenvalue and is the corresponding eigenvector. The 

eigenvalue/eigenvector assignment problem is to determine the number of eigen- 
values in equation (4) that can be arbitrarily assigned and to determine the 
freedom available in the selection of the associated eigenvectors. 


In order to see what freedom exists in the choice of eigenvectors, write 
equation (4) in partitioned form as 


r. 


V.'- 


All 

. Ai2 


■ Bi 


A2I 

. A22_ 

+ 

_B2 

KC 


Svi = XiVi 


(5) 
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I 


where A-|^,Bi ^ Rinxm B-| is nonsingular. Since B is full rank, the non- 

singularity of B-j can be assured, if necessary, by reordering the state vari- 
ables in equations (1). Completing the multiplication of the partitioned 
matrices and some algebraic operations (ref. 8) permits equation (5) to be 
expressed as a set of constraints on the selection of eigenvectors. For clarity 
of presentation, these relations are detailed only for real eigenvalues. Exten- 
sion to complex conjugate pairs in quasi-diagonal form yielding real eigenvector 


pairs is straightforward (ref. 8). 

For real eigenvalues 

[^^^n-m “ = IjG + Xsjz (6) 

[Ai + BiKC]v = Xz (7) 

where X is the eigenvalue, v* = []z' : w'^; v is the eigenvector with 
z s Rfflj and 

S = B2B1-'' (8) 

F = A22 - SA-12 (9) 

G = A21 - SA11 (10) 

Ai = [All : A 12 ] (11) 


Equation (6) represents an underdetermined system of n - m equations in n 
unknowns. Thus m eigenvector entries corresponding to the z-vector can be 
chosen arbitrarily provided X does not coincide with the spectrum of F. 
Examination of equation (7) reveals that at least r eigenvalues and r eigen- 
vectors satisfying equation (6) can be assigned to the system in equations (1) 
by the feedback matrix 

K = Bi-I [z(>")Aj, - AiV(r’)][cv(i’)]"^ (12) 

where is the diagonal matrix of r eigenvalues and and 

have the form = [t-| : t2 J . • . : tj.^ (where ti are vectors). The 

solution to equation (12) is guaranteed provided the eigenvalues/eigenvectors 
are chosen to insure the nonsingularity of [cv(*'^3- It should be noted that 
in the case of state variable feedback (C = In)> all n eigenvalues can be 
assigned to the system provided the modal matrix (matrix of eigenvectors) 

V = [vi : V2 : . . . : vj is nonsingular. An algorithm which constructs such 
a nonsingular V is detailed in reference 8. i^pendix B extends this algo- 
rithm to guarantee the nonsingularity of . 

Notice that by carrying out this analysis on the dual system (B',A',C), 
it can be shown that m eigenvalues can be assigned to the system. This analy- 
sis yields the following well-known result (ref. 9). 
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Lemma 1 ; For system (C,A,B), max(m,r) eigenvalues can be 
assigned using output feedback. 

The analysis so far indicates that only max(m,r) eigenvalues can be 
assigned to the system using output feedback. However, by sacrificing some 
degree of freedom in the selection of the associated eigenvectors it is possible 
to extend the number of eigenvalues that can be assigned to min(n,m + r - 1) 
as is shown in the following section. 


ALGORITHM TO ASSIGN MINIMUM (n,m + r - 1) EIGENVALUES 

The basic approach in the development of this algorithm is^to construct the 
output feedback law in equation (2) as a sum of two feedbacks (K + K). The 
first feedback (K) assigns t eigenvalues, and the second feedback assigns 
additional min(m,n - t) eigenvalues while ensuring the protection of the t 
eigenvalues already assigned. The construction procedure yields a set of suf- 
ficient conditions for assigning min(n,m + r - 1) eigenvalues. These condi- 
tions also help characterize the class of systems which cannot be assigned 
min(n,m + r - 1) eigenvalues. Finally, some design freedom still exists to 
partially assign (r - 1) eigenvectors. 


Step 1 : 


Assign t eigenvalues and corresponding eigenvectors to the system in 
equations ( 1 ) and form the matrix 


N = 


Ni 

’n2 


[VI : V2 ; . . . : vt] 


(13) 


with Ni s and nonsingular. Appendix B details a procedure to construct 
N. Let K be the nonunique feedback (eq. (B6)) corresponding to this assign- 
ment . Then the closed-loop matrix is 


A = A + BKC 


(14) 


Step 2: 

Apply a coordinate transformation 
(C,A,B) (CTi,Ti-‘'aTi,T-i-''b) 


where 


Tl 


Ni . 0 


N2 • In-t 


(15) 
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The transformed system has the form 



'At . Ai 2 ' 




11 

.<*><• 


A 

X + 

• • • 

u 


0 • A 22 


1 

OJ 

<PQ 

1 

> 

y = 

pi : C2]x 





( 16 ) 


where g the diagonal matrix of eigenvalues assigned in step 1. It 

should be noted that C2 is invariant under the transformation T-|. 


Step 3: 


In order to assign additional eigenvalues to the system of equations (16) 
while protecting t eigenvalues (At) already assigned, the second feedback (K) 
is restricted to be of unity rank of the form K = qit' with q ^ R®, ^ R>^, 

and a is chosen so that 


Si' 




[0 : c] 


(17)- 


Now q must be chosen so that 

[a 22 + B2qc] 


( 18 ) 


is assigned min(m,n - t) eigenvalues. 


Since (A22>B2) is controllable, the following result holds: 

A A 

Theorem 1 : The single output subsystem (c,A 22 jB 2) can be 

assigned min(m,n - t) ^eigenvalues if and only if (a) (c,A22) 
is observable and (b) §2 is full rank. 


Theorem 1 follows directly from lemma 1. Further, conditions (a) and (b) 
restrict the admissible set of eigenvalue/eigenvector assignments in step 1 . 

For conceptual convenience, these parametric restrictions are formulated in 
terms of controllability conditions of a fictitious dynamic system in the state 
variable representation. 


After matrix operations in equations (14) and (16) have been performed, 
the submatrix A22 oan be written as 

^22 " [^22 “*■ ^ 2 jj ^12 ® 1^^2 I ( 19 ) 

where N2 = -N2N'j’’^. Transposing equation (19) and writing KC2 = sc, s ^ 
yields 

A 22 = A 22 + c’s^B2 + A^2^2 + c’s’B^N2 (20) 
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Equation (20) can now be looked upon as a system matrix derived by applying 
feedback to the dynamic system 


t = A 22 C + Ai' 2V1 + c'y2 


( 21 ) 


n = B2^ + B-iPi 

with ^ e R^-t, Ti G R“, Pi 
N2 = -N2N1-''. 




P2 e R » Ul = N2C, P2 = s'n, and 


From equations (21) it can be shown that conditions (a) and (b) of theorem 1 
are equivalent to the following conditions: (c) 82 = 82 + N2B-1 and is full 

rank; (d) < j^A22 + controllable. Conditions (c) and (d) clearly 

indicate the restriction on selection of the eigenvector parameters N2 in 
step 1. This yields the following sufficient conditions: 


Theorem 2 : The system (C,A,B) can be assigned 

min(n,m + r - 1) eigenvalues arbitrarily close to the 
desired set if the first t eigenvalues and eigenvectors 
in step 1 are chosen so that 


I N-j is nonsingular 

/V /\/ 


II 82 = B2 + N2B'| and is full rank 
III |[a 22 + ^12^2] is controllable 


"Arbitrarily close" in theorem 2 indicates that slight perturbations in 
eigenvalue specifications may be needed in the following situations: 

(i) Assigned eigenvalues coincide with the spectrum of F (appendix 8). 


(ii) An exact combination of eigenvalue/eigenvector specifications in 
step 1 may not yield a nonsingular N^j (appendix 8) . 

(iii) Coincident spectrum situation similar to situation (i) exists for 
the subsystem eigenvalue assignment of theorem 1 . 


Condition I is required to guarantee the existence of the transformation 
T-j in step 2. This condition can be explicitly included in the synthesis pro- 
cedure as detailed in appendix 8. Condition III is obtained from the property 
that for system (eqs. (21)), the class of feedback from input y-j should be 
restricted so that the controllability of the feedback system with respect to 
the input P2 preserved. Conditions II and III yield nonlinear algebraic 
constraints for the elements of the eigenvector parameter matrix N2 and thus, 
in general, can only be used as test conditions for each assignment in step 1. 
However, example 1 of the section entitled "Numerical Examples" shows how these 
conditions can be explicitly checked. 
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step 4: 


If theorem 2 holds, ^min(m,n - t) eigenvalues can be assigned to the sin- 
gle output system (c,A 22 >B 2 ) by the feedback q derived from lemma 1. 

Step 5: 

The composite feedback law 


u = (K + K)y 

assigns min(n,m + r - 1) eigenvalues to the system. 


( 22 ) 


ASSIGNING n EIGENVALUES 

The development so far has revealed that for systems where n>(m+r-1), 
all system eigenvalues cannot usually be assigned. However, by using equa- 
tions (6) and (7), it is possible to derive conditions for assigning all eigen- 
values of the system as follows. 

If the eigenvalues to be assigned are noncoincident with the spectrum of 
F (if necessary, by a perturbation in specification), then equation (6) can be 
explicitly solved for the w-vector and substituted into equation (7). This sub- 
stitution yields a set of homogeneous equations of the form 

[Hi + B^KLiJzi = 0 (i = 1, 2, . . ., n) (23) 

where Hi and Li are readily derived. Then equation (23) has a nontrivial 
solution if and only if 

rank [Hi + BiKLi] < m for all i (24) 

A set of n nonlinear equations in the m,r parameters of the gain 
matrix K can be derived by setting the appropriate determinant in equation (24) 
equal to zero. However, no general conclusions can be drawn regarding the exis- 
tence of the solution. 

An alternative approach not involving solution of nonlinear equations is to 
assign t eigenvalues as in equation (14) and to attempt to assign the remain- 
ing (n - t) eigenvalues approximately to the subsystem in equation (18) using 
the condition in equation (24). In this case, a set of (n - t) linear equations 
in m unknowns results, and a least-squares solution can be obtained. The 
foregoing discussions assume that r k m if necessary by considering the dual 
system. 
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NUMERICAL EXAMPLES 


Example 1 


Consider the system described in reference 7 where 

To 1 0 o 1 ^ 


A = 


0 

0 


0 

0 


0 

0 


0 0 0 


0 

1 

0 



(25) 


Reduce the system (25) to the special form of equations (1) by ordering state 
variables as (x-| ,x 2 ,X 2 ,xi|) to make Ca_ nonsingular (appendix A) and by apply- 
ing the coordinate transformation (C,A,B) (CTq ,T o~^ATo ,To”^B) where 


10 0 0 


To 


0 111 
0 0 10 


0 0 0 1 


(26) 
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I 


to yield 


A = 


0 

0 

0 

0 


0 

0 

0 

0 


1 

0 

0 

0 




0 

1 

0 

0 


B = 


0 

-1 

1 

0 


C = 


0 
-1 

0 
1 _ 

0 0 0 

1 1 1 


(27) 


V 


From equation (6) it can be shovm that the closed -loop eigenvectors satisfy 


X 0 
0 X 


"2 

ni 


-X (1 - X) 
1 0 


ns 

n4 


( 28 ) 


where the eigenvalue is X and the eigenvector is 

V = (ni n2 ns n4)' 


(29) 


From theorem 2, condition I, n-| ^ 0; therefore, choose n-| = 1 without loss 
of generality. Then, N 2 = (-n 2 -ns -n 4 )'. Condition II is met for all 
choices of N 2 since = 0. Condition III implies that 








■A 


■ 0 

0 

0 ■ 


"r 



-n2 

-n3 

-n4 

f 

1 



1 

0 

0 _ 


1 

( 







(30) 


should be controllable. For the controllability matrix in equation (30) to be 
of full rank, the eigenvector parameters in equation (28) must satisfy 


n2tts - n4 + ns^ + nsn4 f 0 


(31) 


By direct substitution equation (31) is seen to be violated by all admis- 
sible selections in equation (28). Thus, the system cannot be assigned 
(m+r-1)=3 eigenvalues. Indeed, only two eigenvalues can be assigned to 
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this system (ref. 7)- However, the previous analyses (refs. 3 and 7) do not 
precisely lead to this conclusion. 


Example 2 


Consider the system described by (ref. 7) 

To 1 oil 


A 


0 0 1 


0 0 0 


1 


B = 1 



1 1 


c = 


1 0 
0 1 


0 

0 

J 


(32) 


It is required to assign eigenvalues close to -1, -2, and -5. The system (32) 
can be reduced to the form of equations (1) by following the procedure in appen- 
dix A as 


Step 1: Measurement matrix is in the desired form with nonsingular. 

Step 2: Apply coordinate transformation 

(C,A,B) ^ (010,10-^10,10-^6) 


where 


1 0 0 


To = 


0 0 1 


(33) 
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yielding 


A = 


B = 


1 1 

0 1 

0 0 

o' 

-1 
1 




(34) 


C = 


0 1 1 




The system (34) is in the required canonical form of equations (1). The closed- 
loop eigenvector constraints in equations (34) can be derived by identifying the 
respective matrices in equations (8) to (10) as 

I 


Bl = 

_0 

-1 

S = (1 

-1) 


F = 0 


G = (0 -1) 

and equation (6) yields 


Xw = [X 


-(1 + X)] 


zi 

Z2 


(35) 


Equation (35) indicates that z-| and Z 2 can be arbitrarily chosen, provided 
X ^ 0. Now, applying the algorithm to assign (m + r - 1) eigenvalues yields 
the following synthesis sequence. 

Step 1 : 

Assign X = -1. From appendix B, case II, = -^1^, and this implies 

z-j ^ 0 in equation (35). (Condition I, theorem 2.) One acceptable assignment 
is 
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N = (1 1 D’ 


(36) 


and a nonunique feedback gain corresponding to this assignment from equa 
tion (B6) is 



(37) 


The closed-loop system matrix of equation (14) is 



Step 2: 

Transform the system to the canonical form of equations (16) using 



(38) 


(39) 
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From equations (21 ) we can form the dynamic system 


with 



0 0 ' 


T 


'i‘ 



_ ■' 

K * 

J. 

Pi + 

1 

^2 


n = 


( 41 ) 


0 

1 " 

K + 

1 

-1 

i_ 


0 





= (1 

1 )?, 

P 2 


Pi 


J 


Since feedback from Pi does not affect the controllability of system (41) 
with respect to y2> subsystem is pole assignable from condition III of 
theorem 2. Further, the assignment in step 1 satisfies condition II since §2 
is nonsingular. Thus theorem 2 holds. 


Step 3s 


Choose = (-2 1). This choice protects X = -1 assigned in step 1. 

It now remains to choose q so that K = qH' assigns the eigenvalues X2 = -2 
and X3 = -5. 

Step 4: 


Assign X2 = -2, X3 = -5 to the single output subsystem in equation (I8) 

where 


A22 = 



A 


0 

0 


B2 


-1 



0 1 


c = [1 1] J 


(42) 
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with feedback 



Step 5: 

The feedback law of the form of equation (22) is given by 



and assigns the eigenvalues -1, -2, and -5 to the system (eq. (32)). 


(43) 


(44) 


AIRCRAFT LATERAL CONTROL DESIGN 

The advantage of combined control of closed-loop eigenvalues and eigenvec- 
tors using state variable feedback has been investigated in reference 10. The 
utility of the output feedback extensions developed in this report will now be 
illustrated through the design of a lateral controller for an aircraft. 

The linear perturbation model for the lateral motions of an aircraft can 
be modeled as 



(45) 


where x is the state vector of roll rate p, yaw rate 4*) sideslip 3, and 
bank angle 4, respectively. The control vector of aileron 5^ rudder 

angular deflections is u. Roll rate p, yaw rate ijj, and lateral accelera- 
tion ay constitute the output vector y. All angles are in degrees, rates in 
deg/sec, and acceleration in m/sec^ (ft/sec^). 

The respective matrices in equations (45) for a fighter aircraft at an alti- 
tude of 6096 m (20 000 ft), a Mach number of 0.67, and an angle of attack of 
3.45° are given by 
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'-3.79 

0.04 

-52 


1 

-0.14 

-0.36 

4.24 

0 

0.06 

-1 


-0.27 

0 

.05 

1 

0.06 

0 

0 J 

' 25 

9.83" 





1 .42 

-4.2 





0.01 

0.05 





. 0 

0 





1 

0 


0 

o' 


0 

1 


0 

0 


_-0.13 

-0.06 

-3.42 

0 


■ 0 

0 

1 




0 

0 





1 .03 

-0.27 






where the elements of the matrices are approximated to two significant digits. 


Since the output vector in equations (45) is derived as a linear combina- 
tion of both state variables and control inputs, the closed-loop system after 
applying feedback u = K*y takes the form 


X = (A + BKC)x + BPu 


y = (C + DKC)x + DPu 


(46) 


where P = [im - and K = PK* is the equivalent output feedback matrix 

obtained by setting 6=0. Thus, the algorithms developed earlier for systems 
with D = 0 are applicable to systems of the form of equations (45), provided 
P exists. 

Then the feedback gain K* is computed by using the relation 


K* = k[Iij + Dk]-'' 


(47) 
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Direct matrix manipulations show that the inverse in equation (47) exists if P 
exists. 

The handling qualities specifications (ref. 11) imply that the lateral air- 
craft dynamics should be composed of two weakly coupled subsystems. Roll rate 
and bank angle constitute the first subsystem and display predominantly the roll 
subsidence and spiral modes. The second subsystem is characterized by a well- 
damped Dutch roll mode defining the yaw rate and sideslip motions. These speci- 
fications can now be formulated as an eigenvalue/eigenvector assignment problem. 

Table I summarizes the modal characteristics of the free aircraft. From 
the table it is seen that the Dutch roll mode is very lightly damped and appears 
dominantly in the response of the roll variables p and (fi as evidenced by 
dominant entries in the corresponding eigenvector pair. Thus, the eigenvalue/ 
eigenvector modification requires that the closed-loop system have the modes 
and mode-variable associations of table II. (In tables I to III, j = V^.) 


TABLE I.- MODAL CHARACTERISTICS OF FREE AIRCRAFT 



Eigenvalue of - 


Eigenvector 

components 

-3.70 

-0.35 + J2.66 

-0.03 . 


(Roll subsidence) 

(Dutch roll) 

(Spiral) 

P 

-0.964 

-0.403 0.829 

-0.032 


-.041 

-.096 -.131 

.044 

B 

-.002 

.069 -.034 

.002 

4> 

.261 

.324 .112 

.998 


TABLE II.- DESIRED MODAL SPECIFICATIONS 


Mode 

Eigenvalue 

Dominant response variable 

Roll subsidence 

-6 

Roll rate 

Dutch roll 

-1 .0 + jO.2 

Yaw rate and sideslip 

Spiral 

-0.01 

Bank angle 
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The output feedback analysis in the main text shows that all system eigen- 
values can be assigned since n = m + r - 1 and only two eigenvectors (t) can 
be assigned with at most two (m) entries in each vector arbitrarily chosen. The 
eigenvector freedom available was used to control the structure of the eigenvec- 
tor pair corresponding to the Dutch roll mode to effect the desired yaw rate 
and sideslip dominance. The modal coupling matrices (eq. (B4)) aid in 

the selection of the Dutch roll mode to yield the appropriate eigenvector forms. 
The four eigenvector entries that were freely chosen corresponded to the roll 
rate and yaw rate components of the real eigenvector pair associated with the 
Dutch roll mode. Since gain magnitude constraints cannot be explicitly included 
into the synthesis algorithm, the design parameters have to be iteratively modi- 
fied to meet gain limit requirements. After some design iterations, a compro- 
mise design yielded the modal characteristics summarized in table III. For 
example, it was noted that the Dutch roll mode damping could not be reduced (to 
improve and B responses) without violating feedback gain limits which were 
set at unity for this analysis. The design yielded a feedback gain matrix K* 
(eq. (47)) as 


K 




-0.19 “0.2 -0.66 

-0.16 0.19 -0.6 


TABLE III.- MODAL CHARACTERISTICS OF FEEDBACK AUGMENTED AIRCRAFT 


Eigenvector 

components 

Eigenvalue of - 

-6 

(Roll subsidence) 

-1.0 + jO.20 
(Dutch roll) 

-0.011 

(Spiral) 

P 

0.986 

0.01 0.02 

-0.013 


.009 

.35 .53 

-.031 


o 

o 

1 

.25 .72 

-.082 


-.165 

-.05 .05 

-.996 


Table III illustrates that the closed-loop eigenvectors have approached 
the desired mode-decoupled structure. In particular, the desired modification 
achieved in the eigenvector pair corresponding to the Dutch roll mode should be 
noted. The improvement in transient response characteristics using the feedback 
controller is illustrated in figure 1 . The response curves demonstrate that the 
cross coupling between the roll axis (p,4>) and yaw axis (^,B) has been signifi- 
cantly reduced. 
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CONCLUDING REMARKS 


New sufficient conditions to assign minimum (n,m + r - 1) eigenvalues by 
means of output feedback have been derived. In general, in addition to the 
assignment of minimum (n,m + r - 1) eigenvalues, (r - *1) eigenvectors can be 
partially assigned with at most m entries in each vector arbitrarily chosen. 
The utility of assigning both eigenvalues and eigenvectors for response modifi- 
cation is illustrated by designing a feedback controller for the lateral dynam- 
ics of an aircraft. The synthesis algorithm is computationally simple and 
involves only the solution of a linear system of equations. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
December 19, 1977 
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APPENDIX A 


REDUCTION OF (C,A,B) TO CANONICAL FORM 
Consider the given system of the form 



(AD 


Step 1 : 

Reorder the state variables (if necessary) so that the measurement matrix C 
is of the form 

C □ [Ca : Cb] (A2) 

where Cg ^ and nonsingular. It should be noted that the dimensions 

of the partitioned matrices in equation (A2) are different from those in 
equations (3). 

Step 2: 

Apply a coordinate transformation 

X = Tqx (A3) 


where 


Ca~^ . -Ca ^ ( Cb - Cg) 


To = 


(A4) 


and Cg e Rr><n-r and Cg = 


Tn-r 

0 


with c as defined in equations (3). The 


transformed system 

X = (To-‘'aTo)x + To-''Bu| 
y = CTqx 

has the desired form of equations ( D 


(A5) 
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APPENDIX B 


EIGENVALUE/EIGENVECTOR SELECTION PROCEDURE 

For completeness of presentation, an existing algorithm (ref. 8) developed 
for eigenvalue /eigenvector assignment using state variable feedback is detailed 
in this appendix. Extensions are made to adopt the algorithm for the output 
feedback cases discussed in the main text. 


SPECTRAL SYNTHESIS ALGORITHM 


A direct way of constructing a nonsingular modal matrix is to generate the 
eigenvectors which satisfy equation (6) sequentially and insure that they do 
not lie in the eigenspace generated by the vectors already synthesized. The 
algorithm presented accomplishes this construction while constantly testing to 
insure that the set of eigenvectors is a linearly independent set to a degree 
determined by a numerical tolerance parameter set by the designer. For clarity 
of presentation the algorithm is detailed for real eigenvalue assignments. The 
following notations are used throughout the algorithm presentation: 


(1) ej is an n-vector with jth entry equal to 1 and all other entries 
equal to zero . 

(2) vj^ o : wj^J is the kth eigenvector, where zj^. is designer 
specified. 

(3) where i ^ 0, is defined in 

(7) which follows. 


(4) is the jth entry of and for convenience of nota- 

tion is denoted 


(5) is the set of indices (1,2, . . . , n) . 

(6) an elementary upper triangular matrix of order n and index 1, 
takes the form 


with ejO = 0 


= In - ®®j 

and 0 is chosen so that 
M(k)vij(k-1) = a^ej, j e {a(^)^ 


(B1) 


(B2) 


where A^^) is a subset of A^^) (defined in (5)) containing the indices not 
already used in the construction of the matrices . . ., 

can be constructed if and only if ^ 0 (ref. 8) . 




(B3) 


23 



APPENDIX B 


The algorithm is now accomplished by completing steps 1 to 4 for 
k n 1, 2, . . n and by computing the feedback gains from equation (12) 
using C = Ijj. 

Step 1 : 

For A = Xjj, compete = det j^AjjI n-m ■ 

(a) if = 0, perturb Ai^ to j^Ajj + eA^J and repeat the calculation of 

6k . 

(b) if 6k # 0, proceed to step 2. 

Step 2; 

Compute 

= [Xkin-m - ?]■■'[« + XkS] (B4) 

Step 3: 

For some j g {a(w>, 

(a) Compute where 

is the jth row of the transformation matrix q( 1^- 1 ) j g -j ( ^ ) G R® and 
h.(k-l) G Rlx(n-m). 

vl 

(b) Compute 

Gk = [sj^^^J'zk (B5) 

(i) if Ok ^ 0, compute w = D^^^Zk and (from eq. (B2)) and 

return to step 1 for the next value of k. 

(ii) if Ck = 0, select another j g and return to step 3(a). 

(iii) if CTk = 0 for all j e , go to step 4. 

Step 4: 

For some j g {A(k)> , 

(a) if ^ 0, perturb Zk to (zk + szk^ make Ok 0, compute 

Wk and and return to step 1 for the next value of k. 

(b) if = 0> select another j g and return to step 4(a). 

(c) if ° 0 j ^ » perturb Ak to (Ak + eAk) and 

return to step 1 . 
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APPENDIX B 


The following observations can be made regarding the algorithm as outlined. 

(i) Extension of the algorithm to assign complex eigenvalues is straight- 
forward (ref. 8). 

(ii) If some of the eigenvalues of A coincide with those of F, either 
a perturbation can be made in the eigenvalue specification as in step 1 , or 
special eigenvector structures can be derived by noting that equation (6) has 
a solution for Zj^ = 0 (ref. 8). 

(iii) The iterative procedure in part ii of step 3(b) attempts to meet 
exact eigenvalue/eigenvector specifications. In step 4(b) an attempt is made 
to meet exact specifications with slightly relaxed specifications. 

The test in step 4(c) indicates that the Xy. specification implied that the 
corresponding v^^ lies in the eigenspace already generated. The test thus 
demands an eigenvalue perturbation. 

(iv) The Oy^ provide a good measure of the linear independence between 
eigenvectors, if all vectors are normalized to a standard basis. Since the 
determinant of the modal matrix V is given by the product of the 
(k □ 1, 2, . . . , n) , the numerical ill conditioning of V for inversion can 
be effectively controlled by setting a tolerance on the Oyr, 


OUTPUT FEEDBACK EXTENSIONS 

Since the distribution matrix C is already in the special canonical form, 
generating linearly independent vectors to guarantee existence of feedback 
matrices in equations (12) and (14) is simply achieved by restricting the admis- 
sible pivotal indices (eq. (B2)) in the spectral synthesis algorithm as 

follows . 


Case I: Computing K, Equation (12) 

In equation (5), let A^^^^ be the set of r column indices correspond- 
ing to the r linearly independent columns of Cd^). Then apply the spec- 
tral synthesis algorithm using A^(^) instead of A^^) (notation (5)) for 
k □ 1, 2, . . ., r. This application guarantees the invertibility of CV^^) 
in equation (12). 


Case II: Computing K, Equation (14) 

In equation (5), let A^^^^ be the set of t column indices not contain- 
ing the vector c (eqs. (3)). Then apply the spectral synthesis algorithm using 
At^^^ instead of A^ ^ ) (notation (5)) for k=1,2, . , .,t. This applica- 
tion guarantees the nonsingularity of N'l (eq. (13)). Further, since in this 
case CN is rank t, the feedback matrix K (eq. (14)), computed using the 
relation 


KCN = Bi-I [z(t)At Ain] 


(B6) 


APPENDIX B 


to assign t eigenvalues, is not unique. However, a minimum norm least-squares 
solution for K can be computed using the pseudo inverse of CN. 
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(a) p(0) = 10 deg/sec. 

Figure 1.- Comparison of free aircraft and augmented aircraft response to roll 
rate and sideslip step disturbances. A indicates free aircraft response 
and B corresponds to augmented aircraft response. 




0) 

ci 

u 

I 


o 

d> 

be 

0) 

XS 


25 

0 


-25 





(b) 3(0) = 1QO. 

Figure 1 . - Concluded . 
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